DESY 10-083 
June 2010 



ISSN 0418-9833 



Measurement of Charm and Beauty Jets in Deep Inelastic 

Scattering at HERA 



HI Collaboration 



Abstract 

Measurements of cross sections for events with charm and beauty jets in deep inelastic 
scattering at HERA are presented. Events with jets of transverse energy Ed? > 6 GeV and 
pseudorapidity —1.0 < rf ct < 1.5 in the laboratory frame are selected in the kinematic 
region of photon virtuality Q 2 > 6 GeV 2 and inelasticity variable 0.07 < y < 0.625. 
Measurements are also made requiring a jet in the Breit frame with E^ ct > 6 GeV. The 
data were collected with the HI detector in the years 2006 and 2007 corresponding to an 
integrated luminosity of 189 pb _1 . The numbers of charm and beauty jets are determined 
using variables reconstructed using the HI vertex detector with which the impact parame- 
ters of the tracks to the primary vertex and the position of secondary vertices are measured. 
The measurements are compared with QCD predictions and with previous measurements 
where heavy flavours are identified using muons. 
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1 Introduction 



The production of heavy flavour quarks in deep inelastic scattering (DIS) at the HERA electron- 
proton collider is of particular interest for testing calculations in the framework of perturbative 
quantum chromodynamics (QCD). The process has the special feature of involving two hard 
scales: the square root of the photon virtuality Q and the heavy quark mass m. In the case 
of jet production the transverse energy Et of the jet provides a further hard scale. In leading 
order (LO) QCD, the photon-gluon fusion (PGF) processes ep — > eccX and ep — > ebbX axe the 
dominant production mechanisms for charm (c) and beauty (b) quarks respectively. 

The inclusive c and b quark cross sections and the derived structure functions have been 
measured in DIS at HERA using the 'inclusive lifetime' technique [1,2] and found to be well 
described by next to leading order (NLO) QCD. Measurements of the charm cross section using 
the technique of D meson tagging have also been made [3,4] and are found to be in good agree- 
ment with those using inclusive lifetime information. Measurements of the total charm and 
beauty cross sections have been made by identifying their decays to muons [5]. In the charm 
case these measurements show good agreement with the data extracted using the inclusive life- 
time technique, but are somewhat larger in the case of beauty. 

Measurements of beauty quark production using muon tagging have also been made for DIS 
events containing a high E T jet in either the Breit frame [6,7] or in the laboratory frame [8]. As 
in the muon inclusive case [5] the results were found to be somewhat higher than NLO QCD 
predictions, in particular at low values of Q 2 . In photoproduction, measurements of beauty have 
been made using various lepton tagging techniques and have been found to be either somewhat 
higher than [9] or in agreement with [ 1 0] NLO QCD. A measurement in the Breit frame of the 
production of D* mesons in association with high E T dijets [11] was found to be in agreement 
with NLO QCD predictions within the statistics of the measurement. A measurement of c and 
b jets in photoproduction has been made [ ], which uses a similar method to distinguish heavy 
flavour jets as in the present analysis. The results were found to be in good agreement with 
NLO QCD. 

This paper reports on measurements of the cross sections for events with a c or b jet in DIS at 
HERA. The analysis uses an inclusive lifetime technique following a similar procedure as used 
in [1] to distinguish the jets that contain c or b flavoured hadrons from those containing light 
flavoured hadrons only. The data are analysed in the laboratory frame of reference to match 
the acceptance of the HI detector and a heavy flavour jet with the highest transverse energy 
Ejt* > 6 GeV is required. The measurements in the laboratory frame are compared with b 
quark production measurements obtained from muon tagging [ ]. The analysis is extended to 
the Breit frame of reference requiring a jet with transverse energy of E^ ct > 6 GeV. The results 
are also compared with b quark measurements obtained from muon tagging [7]. The cross 
section measurements in both frames of reference are compared with an NLO QCD program 
employing mass factorisation [13]. 

The data for this analysis were recorded in the years 2006 and 2007 with integrated lumi- 
nosities of 135 pb" 1 taken in e + p mode and 54 pb -1 taken in e~p mode . The ep centre of 
mass energy is y/s = 319 GeV, with a proton beam energy of 920 GeV and electron 1 beam 

'In this paper the term 'electron' also denotes 'positron' unless explicitly stated. 
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energy of 27.6 GeV. The measurements are made for the kinematic region of photon virtuality 
Q 2 > 6 GeV 2 and inelasticity variable 0.07 < y < 0.625. 

Jets containing heavy flavoured hadrons are distinguished from those containing only light 
flavours using variables reconstructed using the HI vertex detector. The most important of these 
inputs are the transverse displacement of tracks from the primary vertex and the reconstructed 
position of a secondary vertex in the transverse plane. Hadrons from heavy quark decays typ- 
ically have longer lifetimes than light hadrons and thus produce tracks that have a significant 
displacement from the primary vertex. For jets with three or more tracks in the vertex detector 
the reconstructed variables are used as input to a neural network to discriminate beauty from 
charm jets. 



2 Monte Carlo Simulation 

Monte Carlo simulations are used to correct for the effects of the finite detector resolution, 
acceptance and efficiency. The Monte Carlo program RAPGAP [14] is used to generate DIS 
events for the processes ep — > ebbX, ep — > eccX and ep — > eqX where q is a light quark of 
flavour u, d or s. RAPGAP combines 0(a s ) matrix elements with higher order QCD effects 
modelled by parton showers. The heavy flavour event samples are generated according to the 
massive photon gluon fusion (PGF) matrix element [15] with the mass of the c and b quarks set 
to m c = 1.5 GeV and mj = 4.75 GeV, respectively. The DIS cross section is calculated using 
the leading order 3-flavour parton density function (PDF) set MRST2004F3LO [16]. 

The partonic system for the generated events is fragmented according to the Lund string 
model [17] implemented within the PYTHIA program [18]. The c and b quarks are hadronised 
according to the Bowler fragmentation function [19] using the parameters a = 0.4 GeV" 2 , 
b = 1.03 GeV^ 2 and rq = 1. The HERACLES program [ : ] calculates single photon radiative 
emissions off the lepton line, virtual and electroweak corrections. 

PYTHIA is used to simulate the background contribution from photoproduction jp — > X. 
The assumed heavy flavour cross sections are in agreement with the measurements made by 
HI [12]. 

The samples of events generated for the uds, c, and b processes are passed through a detailed 
simulation of the detector response based on the GEANT3 program [21], and through the same 
reconstruction software as is used for the data. 



3 QCD models 

The jet cross section data in this paper are compared with two approaches within QCD: 

Firstly, the data are compared with the predictions of Monte Carlo programs based on lead- 
ing order matrix elements with the effect of higher orders modelled by initial and final state 
parton showers. The predictions from the RAPGAP Monte Carlo program are calculated with 
the same settings as described in section 2. The renormalisation and factorisation scales are set 
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to // r = /// — Q. The Monte Carlo program CASCADE [22] is also used to produce predictions 
for the b and c jet cross sections. CASCADE is based on the CCFM [23] evolution equation 
and uses off shell matrix elements convoluted with kr unintegrated proton parton distributions. 
The CA SCADE predicti ons use the AO PDF set with m c = 1.5 GeV and m b = 4.75 GeV, and 
/j, r = \/Q 2 + Pt + 4:m 2 , where px is the transverse momentum of the heavy quark in the virtual 
photon-proton centre of mass frame. Due to the fact that the predictions are based on leading 
order matrix elements the uncertainty on the normalisation of the cross sections is large, and is 
not quantified here. 

Secondly, the data are compared with the predictions of the NLO QCD program 
HVQDIS [13]. The program is based on the fixed flavour numbering scheme (FFNS) which 
uses the massive PGF 0{a 2 s ) matrix element [24] and provides weighted events with two or 
three outgoing partons, i.e. a heavy quark pair and possibly an additional light parton. The cal- 
culations are made using the same settings for the choice of the quark masses as for the Monte 
Carlo programs above: m c = 1.5 GeV, m& = 4.75 GeV. At NLO the predictions of QCD de- 
pend on the choice of the scales // r and ///.To investigate the dependence of the predictions on 
the scales two example choices are made. Firstly, the scale /i r = /// = (Q 2 + p\ + m 2 )/2, 
where p T is the transverse momentum of the heavy quark with the highest value of p T in the 
virtual photon-parton centre of mass frame, is used. This choice of scale is motivated by the 
comparison of NLO QCD with recent measurements of inclusive jet data by HI [25]. Secondly, 
the scale // r = /// = \/Q 2 + 4m 2 is selected. This scale has been used in the comparison 
of HVQDIS with HI inclusive and dijet D* DIS data [4, 11]. Since HVQDIS provides cross 
sections at the parton level, corrections to the hadron level are needed in order to compare to 
the data. These corrections are calculated using the RAPGAP Monte Carlo event generator. In 
each kinematic bin of the measurement, the ratio Chad of the RAPGAP hadron level to parton 
level cross sections is calculated and applied as a correction factor to the NLO calculation. The 
hadron level corrections generally amount to a change in the prediction by < 6 % for charm and 
< 15% for beauty. 

In QCD fits to global hard- scattering data the parton density functions are usually extracted 
using the general mass variable flavour number scheme (GM VFNS) [26-32] for heavy quarks. 
This scheme, which interpolates from the massive approach at low scale values to a 'massless' 
approach at high scale values, provides a theoretically accurate description of heavy flavour 
production. Recently, a set of PDFs [ 13] compatible with the FFNS were generated, using the 
standard GM VFNS PDFs [34] to facilitate comparison of the heavy flavour final state data with 
up-to-date PDFs. 

Predictions are made using three different sets of PDFs: the MSTW08FF3 [33] set extracted 
using the GM VFNS but evolved using the FFNS in order to be compatible with HVQDIS; the 
CTEQ5F3 [35] set extracted using the FFNS; and with the CTEQ6.6 [28] set extracted using the 
GM VFNS. The CTEQ6.6 PDF set uses a variable flavour definition of the running coupling a s 
which is different to the fixed flavour definition assumed in HVQDIS. However, the inaccuracy 
introduced by this incompatibility is likely to be compensated by using an up-to-date PDF 
set [36]. 

As an estimate of the uncertainty on each of the NLO QCD predictions the scales // r and 
/// are varied simultaneously by factors of 0.5 and 2, m c is changed by ±0.2 GeV and 
is changed by ±0.25 GeV. The uncertainty from fragmentation is estimated by replacing the 
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Bowler [19] function by the symmetric function in the Lund model [37], corresponding to 
r Q = 0- 

4 HI Detector 

Only a short description of the HI detector is given here; a more complete description may be 
found elsewhere [38, 39]. A right-handed coordinate system is employed at HI, with its origin 
at the nominal interaction vertex, that has its Z-axis pointing in the proton beam, or forward, 
direction and X (Y) pointing in the horizontal (vertical) direction. The pseudorapidity is related 
to the polar angle 9 by 77 = — lntan(0/2). 

Charged particles are measured in the central tracking detector (CTD). This device consists 
of two cylindrical drift chambers interspersed with orthogonal chambers to improve the in- 
coordinate reconstruction and multi-wire proportional chambers mainly used for triggering. 
The CTD is operated in a uniform solenoidal 1.16 T magnetic field, enabling the momentum 
measurement of charged particles over the polar angular range 2 20° < 9 < 160°. 

The CTD tracks are linked to hits in the vertex detector, the central silicon tracker (CST) [40], 
to provide precise spatial track reconstruction. The CST consists of two layers of double-sided 
silicon strip detectors surrounding the beam pipe, covering an angular range of 30° < 9 < 150° 
for tracks passing through both layers. The information on the Z-coordinate of the CST tracks 
is not used in the analysis presented in this paper. For CTD tracks with CST hits in both layers 
the transverse distance of closest approach (DCA) to the nominal vertex in X-Y, averaged over 
the azimuthal angle, is measured to have a resolution of 43 fxm © 51 /im/(P r [GeV]), where 
the first term represents the intrinsic resolution (including alignment uncertainty) and the sec- 
ond term is the contribution from multiple scattering in the beam pipe and the CST; P T is the 
transverse momentum of the particle. The efficiency for linking hits in both layers of the CST 
to a CTD track is around 84%. The efficiency for finding tracks in the CTD is greater than 95%. 

The track detectors are surrounded in the forward and central directions (4° < 9 < 155°) by 
a finely grained liquid argon calorimeter (LAr) and in the backward region (153° < 9 < 178°) 
by a lead- scintillating fibre calorimeter (SPACAL) with electromagnetic and hadronic sections. 
These calorimeters provide energy and angular reconstruction for final state particles from the 
hadronic system and are also used in this analysis to measure and identify the scattered electron. 

Electromagnetic calorimeters situated downstream in the electron beam direction allow de- 
tection of photons and electrons scattered at very low Q 2 . The luminosity is measured with 
these calorimeters from the rate of photons produced in the Bethe-Heitler process ep — > ep"f. 

5 Experimental Method 
5.1 DIS Event Selection 

The events are triggered by a compact, isolated electromagnetic cluster in either the LAr or 
SPACAL calorimeters in combination with a loose track requirement such that the overall trig- 

2 The angular coverage of each detector component is given for the interaction vertex in its nominal position i.e. 
the position of the centre of the detector. 
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ger efficiency is almost 100%. The electromagnetic cluster with the highest transverse energy, 
which also passes stricter offline criteria is taken as the scattered electron. The Z-position of the 
interaction vertex, reconstructed by one or more charged tracks in the tracking detectors, must 
be within ±20 cm of the centre of the detector to match the acceptance of the CST 

Photoproduction events and DIS events with a hard photon radiated from the initial state 
electron are suppressed by requiring ~ Pz,i) > 35 GeV. Here, E,- L and pz,% denote the 

energy and longitudinal momentum components of a particle and the sum is over all final state 
particles including the scattered electron and the hadronic final state (HFS). The HFS particles 
are reconstructed using a combination of tracks and calorimeter deposits in an energy flow 
algorithm that avoids double counting [41]. 

The event kinematics, Q 2 and y, are reconstructed with the 'eE' method [42], which uses 
the scattered electron and the HFS. In order to have good acceptance for the scattered electron 
in the calorimeters the events are selected in the range Q 2 > 6 GeV 2 . The analysis is restricted 
to 0.07 < y < 0.625 in order to ensure there is a high probability of at least one jet within the 
acceptance of the CST and to reduce the photoproduction background. The position of the beam 
interaction region in X and Y (beam spot) is derived from tracks with CST hits and updated 
regularly to account for drifts during beam storage. 

5.2 Jet Reconstruction 

Jets are reconstructed using the inclusive longitudinally invariant k T algorithm with the mass- 
less Pt recombination scheme and the distance parameter i? = 1 in the 77 — plane [43]. 
The algorithm is first run in the laboratory frame using all reconstructed HFS particles and 
the resultant jets are required to have transverse energy E? > 1.5 GeV, in the angular range 
— 1.0 < r] ]Ct < 1.5. The r] range is asymmetric since the y range chosen means few jets have 
7] < —1.0. This cut also means that the jets are not near the boundary between the LAr and 
SPACAL calorimeters. Jets are reconstructed from the Monte Carlo simulation using an identi- 
cal procedure to that of the data. 

The Monte Carlo simulation is also used to define hadron and parton level jets before they 
are processed by the simulation of the detector response. Hadron level jets are defined by 
running the same jet algorithm as for reconstructed jets using all final state particles, including 
neutrinos, but excluding the scattered electron. A Monte Carlo jet at the reconstructed or hadron 
level is defined as a l b jet' if there is at least one b hadron within a cone of radius 1 about the jet 
axis in the r) — (j) plane. A jet is defined as a 'c jet' if there is at least one c hadron within the 
same cone and that c hadron does not arise from the decay of a b hadron. Jets that have not been 
classified as c or b jets are called 'light jets'. Parton level jets are defined for the Monte Carlo 
samples and for the NLO calculation by running the same jet algorithm on final state partons. 
A parton level jet is defined as a b jet if there is at least one b quark within a cone of radius 1 
about the jet axis in the 7] — (f) plane. A parton level jet is defined as a c jet if there is at least one 
c quark and no b quark within the cone. 

In order to compare with perturbative calculations a good correlation between the parton 
level and hadron level jets is necessary. A jet with high transverse energy is required in either 
the laboratory frame of reference E 3 ^ > 6 GeV or in the Breit frame E^ ct > 6 GeV. For the 
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analysis in the laboratory frame the cross section is measured as a function of E^ , ?7 jct , Q 2 , 
the number of jets iVj Ct with E^ > 6 GeV and also for the integrated sample. For the analysis 
in the Breit frame the flavour of the jet is defined in the laboratory, as described above, for jets 
in the range E^ > 1.5 GeV and —1.0 < 77 jct < 1.5. For events satisfying this condition all 
the final state particles are then boosted to the Breit frame using the four vector of the scattered 
electron and the value of Bjorken x obtained from x = Q 2 /sy. The jet finding algorithm is rerun 
on the boosted particles. The jets in the Breit frame are required to have a transverse energy 
Ej! et > 6 GeV and to have a pseudorapidity, when boosted back to the laboratory frame, in the 
range —1.0 < rf ct < 1.5. The cross section is measured as a function of E^ ct and Q 2 for the 
selected events. The data, measured as a function of Q 2 , in both reference frames are compared 
with b jet data obtained from muon tagging, after correcting those results for the muon phase 
space and other, smaller, differences between the kinematic ranges of the measurements. 



5.3 Jet Flavour Separation 

In the analysis presented in this paper the flavour of the event is defined as the flavour of the jet 
with the highest E^ in the laboratory. Therefore, the measured cross sections are proportional 
to the number of events with a jet rather than the number of jets in an event. 

The separation of b, c and light jets is only briefly described here. The procedure closely 
follows that described in [1]. The separation is performed using the properties of those tracks 
which are within a cone of radius 1 from the jet axis in the rj — <f) plane. The tracks are recon- 
structed in the CTD and must have at least 2 CST hits and have transverse momentum greater 
than 0.3 GeV. The impact parameter 5 of a track is the transverse DCA of the track to the beam 
spot point. Tracks with 5 > 0.1 cm are rejected to suppress contributions from the decays of 
long-lived strange particles. 

The number of tracks in the jet after these selections is called V trac k. The track significance 
S is defined as S = S/a(5), where a(S) is the uncertainty on 5. If the angle a between the 
azimuthal angle of the jet </>j e t and the line joining the primary vertex to the point of DCA is 
less than 90°, the significance is defined as positive [1]. It is defined as negative otherwise. The 
significances Si, S 2 and S 3 are defined as the significance of the track with the highest, second 
highest and third highest absolute significance, respectively. The selected tracks are also used 
to reconstruct the position of the secondary vertex. 

The jets are separated into three independent samples. For each sample a different distri- 
bution is used to separate the light, b and c jets. The Si distribution is used for jets where 
-^track = 1 or Si and S 2 have opposite signs. The S 2 distribution is used for the remaining jets 
with iV track = 2 or where S 3 has a different sign to Si and S 2 . Generally S 2 has a better dis- 
crimination between light and heavy flavour jets than Si, since the chance of reconstructing 2 
high significance tracks is small for jets where all the tracks arise at the primary vertex. For jets 
with iVtrack > 3 where Si, S 2 and S 3 all have the same sign an artificial neural network (NN) 
is used to produce a distribution that combines several variables in order to provide an optimal 
discrimination between b and c jets. The inputs to the NN are Si, S 2 , S 3 , the significance of 
the transverse distance between the secondary and primary vertex, the transverse momenta of 
the tracks with the highest and second highest transverse momentum, N tTax ±, and the number 
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of reconstructed tracks at the secondary vertex. The NN is trained using a sample of inclusive 
heavy flavour DIS Monte Carlo events, with b events as 'signal' and c events as 'background', 
as described in [1]. The NN output is signed according to the sign of Si. 

The three distributions that are used in the flavour separation are shown in figure 1. It 
can be seen that the distributions are asymmetric, mainly due to the tracks arising from heavy 
flavour decays. The NN output gives absolute values in the range from about 0.2 to 0.95. The 
light jet distribution is approximately symmetric and peaks towards low absolute values; the 
c and b distributions are asymmetric with more positive than negative entries; the b events are 
peaked towards 1, whereas the c events are peaked towards 0. For the Si, S2 and NN output 
distributions the data are well described by the Monte Carlo simulation and the contribution 
from photoproduction is very small. 

Since the Si, S2 and NN output distributions for light jets are nearly symmetric around zero 
the sensitivity to the modelling of the light jets can be reduced by subtracting the contents of the 
negative bins from the contents of the corresponding positive bins. The subtracted distributions 
are shown in figure 2. The resulting distributions are dominated by c jets, with a b jet fraction 
increasing towards the upper end of the distributions. Overall the light jets contribute only a 
small fraction. 

The fractions of events with c, b and light jets in the data are extracted using a least squares 
simultaneous fit to the subtracted Si, S 2 and NN output distributions (as in figure 2) and the total 
number of events after DIS and jet selection. Only those bins in the significance distributions 
which have at least 25 events before subtraction are considered in the fit, since Gaussian errors 
are assumed. The last fitted bin of the significance distributions, which usually has the lowest 
statistics, is made 3 times as wide as the other bins (see figure 2). 

The uds (light), c and b RAPGAP Monte Carlo simulation samples are used as templates. 
The templates are scaled by factors pi, p c and p&, respectively, to give the best fit. The Monte 
Carlo samples are weighted to the equivalent luminosity of the data sample so that the p scale 
factors are the ratio between the cross sections of the Monte Carlo models and the data. The 
PYTHIA Monte Carlo program is used to estimate photoproduction background and found to 
be 0.8% overall. The contributions of light, c and b jets in photoproduction are fixed to the 
PYTHIA prediction. Only the statistical errors of the data and Monte Carlo simulations are 
considered in the fit. 

The fitted p parameters for the whole kinematic range and for each of the differential distri- 
butions are listed in table 1. The table includes the correlation coefficients of the fit parameters. 
The fitted parameter p c is seen to be anti-correlated with both pi and pb, due to c jets being a 
significant contribution to the total jet cross section. The magnitude of the correlation C\ c is 
greater than Cb c reflecting the fact that the shapes of the Monte Carlo templates for c jets are 
more similar to those for the light jets than those for b jets. Also included in the table is the 
X 2 /n.d.f. for each fit evaluated using statistical errors only. Acceptable values are obtained for 
all fits. 

The fitted p c value for each bin is converted to a c jet cross section using 

n /V MC g en 

Oc = , (1) 
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where A^ Cgen is the number of generated events that pass the DIS kinematic selection of the bin 
and which contain a c jet passing the jet cuts of the bin at the hadron level, C is the integrated 
luminosity of 189 pb" 1 and C ra d is a radiative correction, calculated from the HERACLES 
Monte Carlo program. The number of generated events N^ Cgcn is calculated after normalising 
the luminosity of the Monte Carlo samples to that of the data as described above. The b cross 
sections are evaluated in a corresponding manner. The differential cross sections are obtained 
from the cross sections integrated over the bin interval by dividing by the size of the bin interval, 
and no further bin centre correction is applied. 



6 Systematic Uncertainties 

The following uncertainties are taken into account in order to evaluate the systematic error. 

• The uncertainty in the 5 resolution of the tracks is estimated by varying the resolution 
by an amount that encompasses any difference between the data and the simulation. This 
was achieved by applying an additional Gaussian smearing in the Monte Carlo simulation 
of 200 fim to 5% of randomly selected tracks and 12 /im to the rest. 

• The uncertainty due to the track efficiency uncertainty is estimated by varying the effi- 
ciency of the CTD by ±1% and that of the CST by ±2%. 

• The uncertainties on the various D and B meson lifetimes, decay branching fractions and 
mean charge multiplicities are estimated by varying the input values of the Monte Carlo 
simulation by the errors on the world average measurements. For the branching fractions 
of b quarks to hadrons and the lifetimes of the D and B mesons the central values and 
errors on the world averages are taken from [44]. For the branching fractions of c quarks 
to hadrons the values and uncertainties are taken from the e + e~ average of [45], which are 
consistent with measurements made in DIS at HERA [46]. For the mean charged track 
multiplicities the values and uncertainties for c and b quarks are taken from Marklll [47] 
and LEP/SLD [48] measurements, respectively. 

• The uncertainty on the fragmentation function of the heavy quarks is estimated by reweight- 
ing the events according to the longitudinal string momentum fraction z carried by the 
heavy hadron in the Lund model using weights of (1 =f 0.7) • (1 — z) + z ■ (1 ± 0.7) for 
charm quarks and by (1 =)= 0.5) • (1 — z) + z ■ (1 ± 0.5) for beauty quarks. The variations 
for the charm fragmentation are motivated by encompassing the differences between the 
Monte Carlo simulation and HI D* data [49]. The size of the variations is reduced for 
beauty compared with charm since the fragmentation spectrum is harder. 

• The uncertainty on the QCD model of heavy quark production is estimated by reweight- 
ing the jet transverse momentum and pseudorapidity by /(10 GeV)) ±0,2 and (1 ± 
^et)±o.i5 for charm jets and (E^/(10 GeV)) ±a3 and (1 ±^)±o.s f or beauty jets. These 
values are obtained by comparing these variations with the measured cross sections. 
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• The uncertainty on the asymmetry of the light jet 5 distribution is estimated by repeating 
the fits with the subtracted light jet distributions (figure 2) changed by ±30%. The light 
jet asymmetry was checked to be within this uncertainty by comparing the asymmetry 
of Monte Carlo simulation events to that of the data for K° candidates, in the region 
0.1 < \5\ < 0.5 cm, where the light jet asymmetry is enhanced. 

• The uncertainty on the reconstruction of (pj et is estimated by shifting its value by ±2°. 
The uncertainty was evaluated by comparing the distribution of the difference between 
4>j et and the track azimuthal angle in data and Monte Carlo simulation. 

• The uncertainty arising from the hadronic energy scale is estimated by changing the 
hadronic energy by ±2% for jets in the laboratory and ±4% for the jets in the Breit 
frame. 

• The uncertainty arising from the electron energy scale and polar angle is estimated by 
changing the electron energy by ±1% and the polar angle by ±1 mrad. 

• The uncertainty in the photoproduction background is estimated by varying the expected 
number of events by ±100%. 

• The uncertainty on the luminosity is 4%. 

• The uncertainty on the radiative correction is 2%. 

The above systematic uncertainties are evaluated by making the changes described above to 
the Monte Carlo simulation and repeating the procedure to evaluate the c and b cross sections, 
including the fits. The uncertainties are evaluated separately for each measurement bin and are 
treated as correlated except for the radiative corrections. 

The most important sources of systematic error for the charm jets are the uncertainty on the 
light jet contribution, the uncertainty of the impact parameter resolution and the contribution of 
the uncorrelated errors. For the beauty jets, the systematic uncertainties are considerably larger 
with the main sources of uncertainty being those due to the multiplicity of b quark decays, the 
track efficiency, the hadronic energy scale and the impact parameter resolution. 



7 Results 

The cross sections for c and b jets are presented in the laboratory frame of reference (section 7.1) 
and in the Breit frame (section 7.2). The b jet data are also compared with measurements ob- 
tained from muon tagging (section 7.3). The cross sections for events with c or b jets are shown 
together with theoretical predictions in table 2. The cross section values for all the measure- 
ments are given in table 3 with the contribution of the systematic errors for each measurement 
listed in table 4. 
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7.1 Jet Cross Sections in the Laboratory Frame 

The jet cross sections in the laboratory frame are measured in the kinematic range Q 2 > 6 GeV 2 
and 0.07 < y < 0.625 for the heavy flavour jet with the highest Ej, with > 6 GeV and 
— 1.0 < rf ct < 1.5. The hadron level c and b cross sections with jets are 

3290 ± 50(stat.) ± 260(syst.) pb 

and 

189 ± 9(stat.) ± 42(syst.) pb, 

respectively. Here the first error is statistical and the second is systematic. 

These cross sections are compared in table 2 to the expectations of the Monte Carlo pro- 
grams RAPGAP and CASCADE as well as to the NLO predictions with HVQDIS including 
hadronisation corrections. The NLO predictions are given for three different sets of PDFs and 
two different scale choices, /i = a/ (Q 2 + p\ + m?)/2 and /i = \/Q 2 + 4m 2 . Overall, RAP- 
GAP agrees well with data for both charm and beauty. CASCADE predicts a significantly 
larger c cross section than observed in data while for beauty the discrepancy is much reduced. 
Within uncertainties the NLO predictions agree reasonably well with the data both for charm 
and beauty. In general the NLO expectations for beauty display a smaller dependence on scale 
than for charm. 

Differential c and b jet cross sections are measured as a function of E%\ V et , Q 2 , and the 
number of jets N' sct with E^ > 6 GeV (table 3). The differential c cross sections are shown in 
figure 3 in comparison to Monte Carlo expectations. The RAPGAP model describes all these 
distributions reasonably well in shape and normalisation. CASCADE exceeds the data at low 
E 3 t as well as at low Q 2 but provides a good description at high Q 2 and high E? , respectively. 
The excess of CASCADE is concentrated in the forward rf ct region. As expected from the 
visible cross section given in table 2, CASCADE lies above the data in the N^ ct distribution. 
The model does, however, give a reasonable description of the rf ct distribution after accounting 
for the difference in normalisation. 

The charm jet cross section measurements are shown in figure 4 together with the NLO 
predictions of HVQDIS. In general the NLO expectations describe the data reasonably well in 
all differential distributions although the predictions with the scale \i = a/ (Q 2 + p\ + m 2 )/2 
fall somewhat below the data at low Q 2 , low E^ and in the forward rf ct region. 

In figures 5 and 6 the differential b cross sections are shown as a function of Ej?, rf et , Q 2 
and N :,ct in comparison to Monte Carlo and NLO expectations, respectively. RAPGAP yields a 
good description of all distributions also for beauty. CASCADE overshoots the data at small Q 2 , 
is slightly above the data at small Ej, and shows an excess in the forward if et direction. These 
differences are similar but less significant than for charm. HVQDIS gives a good description of 
the beauty data with little dependence on the choice of scale. 

7.2 Jet Cross Sections in the Breit Frame 

Differential c and b cross sections are also measured for the highest E^ e jet in the Breit frame 
with E^ et > 6 GeV in the kinematic range Q 2 > 6 GeV 2 , 0.07 < y < 0.625 for the heavy 
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flavour jet with the highest E^ in the laboratory satisfying E^ > 1.5 GeV and — 1.0 < ?7 jct < 
1.5. 

The c cross sections are shown as a function of Q 2 and Ej! ct in figure 7. The data are com- 
pared to the expectations from RAPGAP, CASCADE and HVQDIS. It can be seen that RAP- 
GAP provides a good description of the data for both distributions. CASCADE overestimates 
the data in Q 2 and E^ ct . In contrast to the observation in the laboratory frame, the deviation in 
Q 2 is found to be independent of Q 2 here. Nevertheless the shapes of the predictions are sim- 
ilar to those in the data. As for the laboratory frame analysis, HVQDIS with the scale choice 
/i = \/Q 2 + 4m 2 reproduces the data well, while for the scale \x = a/ {Q 2 + Pt + rn 2 )/2 it 
tends to underestimate the c jet data at low values of Q 2 and E^ et . 

The differential b cross sections are shown as a function of Q 2 and Eji ct in figure 8 together 
with the Monte Carlo and NLO expectations. As for charm, RAPGAP performs well while 
CASCADE lies systematically above the data. The high rate of Breit frame jets in CASCADE 
both for charm and beauty jet production is related to the transverse momentum distribution 
of the unintegrated gluon density used for the calculations. HVQDIS describes the data well 
showing little dependence on the choice of scales. 

7.3 Comparison with Muon Tagging Measurements 

The b jet cross sections may be compared with b jet measurements obtained from muon tagging 
in the Breit (HI [ ]) and laboratory (ZEUS [8]) frames of reference. The muon measurements 
were made requiring the presence of a muon and a jet in either the laboratory frame, with 
E ] £ > 5 GeV or in the Breit frame with E^ ct > 6 GeV and with a central rapidity requirement 
in the laboratory frame, similar to the present analysis. The measurements were also made in a 
similar y range but start at lower values of Q 2 (Q 2 > 2 GeV 2 ). Therefore, comparison of the 
cross sections with these measurements as a function of E^ or E^ e would require interpolating 
over a large range in Q 2 . However, the b cross sections can be compared as a function of Q 2 for 
the range where the Q 2 binning of the muon measurements overlaps closely with the present 
analysis, namely Q 2 > 10 GeV 2 for the laboratory analysis and Q 2 > 6 GeV 2 for the Breit 
frame analysis. 

The present analysis is repeated with two different sets of Q 2 bins chosen to match the 
HI and ZEUS muon measurements as closely as possible. The cross sections are shown as a 
function of Q 2 for the two sets of bins in figure 9. The HI muon data are corrected by factors 
of about 15 which are obtained using the RAPGAP Monte Carlo. The dominant corrections 
account for the b — > \i branching fraction and for the extrapolation from the phase space of 
the muon measurement, which had restrictions on p^ and rf, to the phase space of the present 
analysis. The ZEUS muon data are corrected to the present phase space by factors of around 
6. These corrections are smaller than in the case of the HI data because the ZEUS data have a 
wider rf and p^ coverage. The corrections also include smaller effects due to the difference in 
the Et range of the jets for the ZEUS laboratory frame analysis, differences in the r\ ranges of 
the jets, the difference in the y ranges, the difference in the jet cross section definitions and jet 
finding algorithms and the fact that the lower edge of the lowest Q 2 bin is Q 2 > 5 GeV 2 for 
the HI muon measurement. An additional uncertainty of around 10% is added to the corrected 
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muon measurements to account for theoretical uncertainties on the extrapolation factors coming 
from uncertainties on the perturbative scales and fragmentation model used. The central values 
of the present data in the Breit frame are found to lie below the adjusted HI muon data at 
high Q 2 . The present data in the laboratory frame are found to lie significantly below the ZEUS 
muon data at low Q 2 , where the difference is a factor 2.1. The comparison suggests a systematic 
difference between the HI inclusive lifetime tagged data and the muon tagged data, particularly 
in comparison with the ZEUS muon tagged data at low Q 2 . 

8 Conclusion 

The cross sections for events with charm and beauty jets have been measured in deep inelastic 
scattering at the HERA electron-proton collider. Measurements are made in the laboratory 
frame for Ej! > 6 GeV and —1.0 < rj < 1.5 for the kinematic region of photon virtuality 
Q 2 > 6 GeV 2 and inelasticity variable 0.07 < y < 0.625. Measurements are also made in the 
Breit frame of reference. The analysis uses the precise spatial information from the HI vertex 
detector to distinguish those jets that contain c and b flavoured hadrons from jets containing 
only light flavoured hadrons. 

The laboratory frame jet data are compared with the Monte Carlo models RAPGAP and 
CASCADE. RAPGAP is generally found to give a good description of the data. CASCADE 
is found to lie above the charm data, especially at low Q 2 and high rf et . After accounting for 
the difference in normalisation CASCADE generally gives a good description of the shape of 
the differential cross section measurements. CASCADE provides a better prediction of the 
beauty cross section normalisation than it does for charm but still tends to overestimate the 
data at low Q 2 and high rf**. The data are also compared with NLO QCD calculations made 
using the HVQDIS program. The beauty data are well described by the calculation. The charm 
expectations are found to depend strongly on the choice of renormalisation and factorisation 
scale. The differential cross sections are described within the experimental and theoretical 
uncertainties with a scale choice of \i = /i r = /i/ = \JQ 2 + 4m 2 . The predictions tend to lie 
below the data at low Q 2 and high rj with a choice of = a/ (Q 2 + p% + m 2 )/2. 

For the measurements of the cross section requiring a jet in the Breit frame with E^ e > 
6 GeV RAPGAP is found to give a good description of the data while CASCADE again over- 
estimates the cross sections. The NLO QCD predictions for charm jets with a scale choice of 
H = \/Q 2 + 4:m 2 are compatible with the data while the predictions with the choice of scale 
/i = a/ (Q 2 + Pj> + m 2 )/2 fail to describe the data at low Q 2 . The b jet data are described by 
NLO QCD for all choices of scale. 

The b jet data are compared with HI and ZEUS data obtained from muon tagging by ad- 
justing that data mainly for the extrapolation of the measured to the full muon phase space and 
for the b — > fi branching fraction. The b jet data from the present analysis are found to lie 
systematically below those obtained from ZEUS at low Q 2 and below the HI muon tagged data 
at high Q 2 . 

The present measurements show that charm and beauty production in deep inelastic scatter- 
ing, adequately described by NLO QCD in the inclusive case, is also described in the presence 
of an additional hard scale provided by a jet. 
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bin 


Q 2 range 
(GeV 2 ) 


7-i(*)jet 

' range 
(GeV) 


^jot range 




Pi 


Pc 


Pb 


X 2 /n.d.f. 




Cib 


Cfcc 


1 


Q 2 > 6 


> 6 


-1.0<77 jGt < 1.5 


> 1 


1.178 ±0.004 


1.040 ±0.015 


0.95 ±0.05 


46.3/49 


-0.95 


0.52 


-0.66 


2 
3 
4 
5 


Q 2 > 6 


6 <£^ ct < 10 
10 <E^ < 16 
16 <E^ < 24 
24 <E^ < 36 


-1.0<r/j ct < 1.5 


> 1 


1.203 ±0.006 
1.142 ±0.009 
1.060 ±0.016 
1.080 ±0.030 


1.028 ±0.020 
1.078 ±0.030 
1.092 ±0.066 
0.713 ±0.158 


0.93 ±0.11 
1.00 ±0.06 
0.87 ±0.09 
0.93 ±0.20 


54.8/46 
34.8/45 
41.0/40 
16.0/30 


-0.95 
-0.95 
-0.95 
-0.96 


0.61 
0.54 
0.51 
0.50 


-0.77 
-0.67 
-0.64 
-0.62 


f. 
U 

7 
8 
9 
10 


D 2 ^ R 




— l.U\/y *\ — u.o 

-0.5<77 jct < 0.0 
0.0<J7 jct < 0.5 
0.5<r/j ct < 1.0 
1.0<77 jct < 1.5 


1 


1.11Z ZC U.U1Z 

1.127 ±0.009 
1.199 ±0.008 
1.213 ±0.009 
1.188 ±0.014 


U.OOO ZC U.UoO 

0.985 ± 0.027 
1.020 ±0.029 
1.172 ±0.033 
1.209 ±0.060 


n qr + n 91 

U.JO ZIZ U.Z1 

0.92 ±0.09 
1.05 ±0.08 
0.87 ±0.08 
0.81 ±0.14 


oa o //in 
ou.y / ^u 

47.5/44 

46.5/45 

37.9/43 

39.8/41 


— u.yo 
-0.95 
-0.94 
-0.95 
-0.97 


u.o i 
0.52 
0.51 
0.50 
0.60 


n 7i 

— U. f 1 

-0.66 
-0.64 
-0.64 
-0.73 


1 1 
1 1 

12 
13 
14 
15 


fi ^n 2 ^is 
18 <Q 2 < 45 
45 <Q 2 < 110 
110 <Q 2 < 316 
316 <Q 2 < 1000 




— i.U\/y *\ i.o 


<S 1 


l.oo 1 in u.U±± 

1.153 ±0.008 
1.091 ±0.007 
1.177 ± 0.011 
1.084 ±0.016 


1 1 7 a n 
i.i i q nz u.uoo 

1.109 ±0.027 

0.986 ±0.028 

0.917 ±0.039 

0.866 ±0.065 


u. i o x u.uy 
1.00 ±0.10 
1.11 ±0.10 
1.06 ±0.10 
0.86 ±0.12 


44.2/43 
44.0/45 
41.2/42 
39.0/39 


— u.yo 
-0.95 
-0.95 
-0.95 
-0.95 


U.OO 

0.54 
0.53 
0.53 
0.50 


— u.o f 
-0.68 
-0.67 
-0.67 
-0.63 


16 
17 


2 > 6 


£^ et > 6 


-1.0<r/j ct < 1.5 


= 1 
= 2 


1.208 ±0.005 
0.918 ±0.013 


1.034 ±0.016 
1.074 ±0.045 


0.98 ±0.06 
0.88 ±0.07 


47.4/49 
42.7/40 


-0.95 
-0.95 


0.53 
0.53 


-0.67 
-0.67 


18 
19 
20 


Q 2 >6 


6 <E^ ct < 10 
10 <£* jot < 16 
16 <£* jct < 24 


-1.0<r/j ct < 1.5 


> 1 


1.068 ±0.009 
0.948 ±0.016 
0.803 ±0.033 


1.111 ±0.029 
0.987 ±0.055 
0.973 ±0.129 


0.86 ±0.08 
1.02 ±0.07 
0.86 ±0.13 


44.7/44 
41.9/41 
33.1/33 


-0.95 
-0.95 
-0.96 


0.56 
0.53 
0.50 


-0.72 
-0.67 
-0.64 


21 
22 
23 
24 
25 


6 <0 2 < 18 
18 <Q 2 < 45 
45 <Q 2 < 110 
110 <Q 2 < 316 
316 <Q 2 < 1000 


E^ et > 6 


— 1 0<77 jct < 1 5 


> 1 


1 276 ± 016 
0.999 ±0.014 
0.864 ±0.014 
0.927 ±0.018 
0.865 ±0.024 


1 248 ± 047 
1.120 ±0.045 
0.949 ± 0.049 
0.945 ± 0.068 
0.775 ±0.102 


76 ± 1 
0.98 ±0.09 
1.04±0.11 
0.95 ±0.12 
0.74 ±0.15 


35 6 /40 
29.7/40 
42.6/39 
38.5/37 
16.2/31 


—0.95 
-0.95 
-0.95 
-0.95 
-0.95 


0.54 
0.54 
0.53 
0.52 
0.48 


—0.69 
-0.69 
-0.68 
-0.67 
-0.62 
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6 <Q 2 < 18 
18 <Q 2 < 100 


£* Jct > 6 


-1.0<r/ jct < 1.5 


> 1 


1.272 ±0.016 
0.944 ±0.010 


1.250 ±0.047 
1.055 ±0.034 


0.77 ± 0.10 
1.01 ±0.07 


36.5/40 
31.9/43 


-0.95 
-0.95 


0.54 
0.54 


-0.68 
-0.69 


28 
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30 


10 <Q 2 < 25 
25 <Q 2 < 100 
100 <Q 2 < 1000 


E* T }Ct > 6 


-1.0<r/j ct < 1.5 


> 1 


1.277 ±0.010 
1.109 ±0.006 
1.131 ±0.009 


1.127 ±0.029 
1.044 ±0.022 
0.919 ±0.031 


0.86 ±0.09 
1.06 ±0.08 
0.96 ±0.08 


42.9/43 
45.9/46 
37.7/44 


-0.95 
-0.95 
-0.95 


0.53 
0.53 
0.51 


-0.67 
-0.68 
-0.65 



Table 1: The fit parameters pi, p c and pb along with their errors, the \ 2 P er degree of freedom and the correlation coefficients. The first 
row lists the results of the fit used to evaluate the integrated cross sections (bin 1). The remaining rows lists the fits used to evaluate the 
differential cross sections for jets in the laboratory frame (bins 2-17 and 28-30) and those requiring at least one jet in the Breit frame (bins 
18-27). 
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r 1 1 
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3290 ± 50 ± 260 
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AO 
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V(Q 2 +P 2 T + m 2 )/2 


MSTW08FF3 


278011° 
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^/(Q 2 +p 2 T + m 2 )/2 


CTEQ6.6 


2780338 


196™ 
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y/(Q 2 +p 2 T + m 2 )/2 


CTEQ5F3 
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^Q 2 + 4m 2 




280011° 
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Table 2: The cross sections for events with c and b jets for the kinematic range Q 2 > 6 GeV 2 , 
0.07 < y < 0.625, > 6 GeV and -1.0 < r/ jct < 1.5. The measured data cross sections 
are shown with their statistical and systematic uncertainties. The data are compared with the 
predictions from the Monte Carlos RAPGAP and CASCADE and with NLO QCD, calculated 
using HVQDIS. The NLO QCD predictions are shown for three sets of parton distribution 
functions and two choices of renormalisation and factorisation scales. The errors are obtained 
by changing the scales by factors of 0.5 and 2, by varying the quark masses and using a different 
model for the fragmentation of the quarks. 
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1.00 


b28 












41.6 


10.7 


23.9 


1.07 


c29 


25 <Q 2 < 100 










1400.1 


2.1 
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1.00 


b29 












77.9 


7.7 


22.7 
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c30 


100 <Q 2 < 1000 










664.0 


3.4 


9.1 


0.99 


b 30 












47.5 


7.9 


19.8 


1.01 



Table 3: The measured charm and beauty cross sections for those events in which the highest 
jet is a charm or beauty jet. Integrated cross sections in each bin are shown. The first 
two rows (bin 1) are the integrated charm and beauty cross sections for the measured phase 
space respectively. The differential cross sections may be formed from the remaining rows by 
dividing by the corresponding bin width. The remaining rows list the cross sections for jets in 
the laboratory frame (bins 2-17 and 28-30) and those requiring at least one jet in the Breit frame 
(bins 18-27). The data is corrected to the hadron level. The table also shows the statistical(5 sta t) 
and systematic error (<5 sys ), together with the hadronic correction Chad that is applied to the NLO 
theory to compare with the data. 
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Table 4: The contributions to the total systematic error. The bin numbering scheme follows that used in table 3. The first column lists 
the uncorrected systematic error. The next 10 columns represent a +la shift for the correlated systematic error contributions from: track 
impact parameter resolution; track efficiency; c fragmentation; b fragmentation; light quark contribution; struck quark angle qua rk; hadronic 
energy scale; photoproduction background; electron energy scale; electron theta; reweighting the jet transverse momentum distribution P^ ct 
and pseudorapidity rf et distribution for c and b events; the c hadron branching fractions and multiplicities; and the b quark decay multiplicity. 
Only those uncertainties where there is an effect of > 1% in any bin are listed separately; the remaining uncertainties are included in the 
uncorrected error. There is an additional contribution to the systematic error due to the uncertainty on the luminosity of 4%. 
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Figure 1 : The significance distribution Si (a), S 2 (b) and the output of the neural network (NN 
Output) (c) for tracks of the highest transverse energy jet in the event. Included in the figure is 
the expectation from the Monte Carlo simulation for uds (light), c and b events. The contribu- 
tions from the various quark flavours in the Monte Carlo simulation are shown after applying 
the scale factors p h p c and pb, as described in the text. The background (BG) contribution from 
a photoproduction Monte Carlo simulation is also shown. 
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Figure 2: The subtracted distributions of S\ (a), S 2 (b) and the neural network output (c) for 
the highest transverse energy jet in the event. Included in the figure is the result from the fit to 
the data of the Monte Carlo simulation distributions of the uds (light), c and b quark flavours to 
obtain the scale factors p t , p c and pb, as described in the text. The background (BG) contribution 
from a photoproduction Monte Carlo simulation is also shown. 
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Figure 3: The differential cross sections for the highest transverse energy charm jet in the 
laboratory frame as a function of , rf et , Q 2 and the number of laboratory frame jets in the 
event N^ ct . The measurements are made for the kinematic range > 6 GeV, — 1 < tf** < 
1.5, Q 2 > 6 GeV 2 and 0.07 < y < 0.625. The inner error bars show the statistical error, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data are 
compared with the predictions from the Monte Carlo models RAPGAP and CASCADE. 
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Figure 4: The differential cross sections for the highest transverse energy charm jet in the 
laboratory frame as a function of E 3 ^ 1 , if et , Q 2 and the number of laboratory frame jets in the 
event iVj ct . The measurements are made for the kinematic range Ej, > 6 GeV, — 1 < 7] }et < 
1.5, Q 2 > 6 GeV 2 and 0.07 < y < 0.625. The inner error bars show the statistical error, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data 
are compared with the predictions from NLO QCD where the bands indicate the theoretical 
uncertainties. 
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Figure 5: The differential cross sections for the highest transverse energy beauty jet in the 
laboratory frame as a function of , rf et , Q 2 and the number of laboratory frame jets in the 
event N^ ct . The measurements are made for the kinematic range > 6 GeV, — 1 < tf** < 
1.5, Q 2 > 6 GeV 2 and 0.07 < y < 0.625. The inner error bars show the statistical error, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data are 
compared with the predictions from the Monte Carlo models RAPGAP and CASCADE. 
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Figure 6: The differential cross sections for the highest transverse energy beauty jet in the 
laboratory frame as a function of , ff et , Q 2 and the number of laboratory frame jets in the 
event iVj ct . The measurements are made for the kinematic range Ej, > 6 GeV, — 1 < 77 jet < 
1.5, Q 2 > 6 GeV 2 and 0.07 < y < 0.625. The inner error bars show the statistical error, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data 
are compared with the predictions from NLO QCD where the bands indicate the theoretical 
uncertainties. 
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Figure 7: The differential cross sections da / dE^ ct and da / dQ 2 for events with a jet in the Breit 
frame, where the jet with the highest transverse energy in the laboratory frame satisfying > 
1.5 GeV and —1 < ?7 jct < 1.5 is a charm jet. The measurements are made for the kinematic 
range Q 2 > 6 GeV 2 and 0.07 < y < 0.625. The inner error bars show the statistical error, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data are 
compared with the predictions from the Monte Carlo models RAPGAP and CASCADE (upper 
plots) and the NLO QCD calculation (lower plots), where the bands indicate the theoretical 
uncertainties. 
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Figure 8: The differential cross sections da / dE^ ct and da / dQ 2 for events with a jet in the Breit 
frame, where the jet with the highest transverse energy in the laboratory frame satisfying > 
1.5 GeV and — 1 < ?7 jct < 1.5 is a beauty jet. The measurements are made for the kinematic 
range Q 2 > 6 GeV 2 and 0.07 < y < 0.625. The inner error bars show the statistical error, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data are 
compared with the predictions from the Monte Carlo models RAPGAP and CASCADE (upper 
plots) and the NLO QCD calculation (lower plots), where the bands indicate the theoretical 
uncertainties. 
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Figure 9: The upper plots show the differential cross section da/dQ 2 for events with a jet 
in the Breit frame with E^ ct > 6 GeV, where the jet with the highest transverse energy in the 
laboratory frame satisfying Ej, > 1.5 GeV and— 1 < if ct < 1.5 is a beauty jet. The lower plots 
show the differential cross section da/dQ 2 for events with a beauty jet in the laboratory frame 
with > 6 GeV and — 1 < rf ct < 1.5. The present measurements are made for the kinematic 
range Q 2 > 6 GeV 2 and 0.07 < y < 0.625. The inner error bars show the statistical error, the 
outer error bars represent the statistical and systematic errors added in quadrature. The data 
are compared with the measurements obtained using muon tagging from HI [7] (upper plots) 
and ZEUS [8] (lower plots) extrapolated to the present phase space and shifted in Q 2 for visual 
clarity. For the muon data the outer error bars show the statistical, systematic and extrapolation 
uncertainties added in quadrature. The data are also compared with the predictions from the 
Monte Carlo models RAPGAP and CASCADE (left) and the NLO QCD calculation (right), 
where the bands indicate the theoretical uncertainties. 
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